WELL-POSEDNESS AND GLOBAL ATTRACTORS FOR LIQUID 
CRYSTALS ON RIEMANNIAN MANIFOLDS 



STEVE SHKOLLER 

Abstract. We study the coupled Navier-Stokes Ginzburg-Landau model of nematic liq- 
uid crystals introduced by F.H. Lin, which is a simplified version of the Ericksen-Leslie 
system. We generalize the model to compact n-dimensional Riemannian manifolds, and 
show that the system comes from a variational principle. We present a new simple proof for 
the local well-posedness of this coupled system without using the higher-order energy law. 
We then prove that this system is globally well-posed and has compact global attractors 
when the dimension of the manifold M is two. Finally, we introduce the Lagrangian aver- 
aged liquid crystal equations, which arise from averaging the Navier-Stokes fluid motion 
over small spatial scales in the variational principle. We show that this averaged system is 
globally well-posed and has compact global attractors even when M is three-dimensional. 
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1. Introduction 

Nematic liquid crystals are well-studied and interesting examples of anisotropic non- 
Newtonian fluids. A liquid crystal is a phase of a material between the solid and liquid 
phases. The solid phase has strong intermolecular forces that keep the molecular position 
and orientation fixed, while in the liquid phase, the molecules neither occupy a specific 
average position nor do they remain in any particular orientation; the liquid crystal phase 
does not have any positional order, but does possess a certain amount of orientational order. 
This phase is described by a velocity field, as well as a director field that describes locally 
the averaged direction or orientation of the constituent molecules. In this paper, we shall 
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analyze the behavior of a certain model of nematic liquid crystals on compact Riemannian 
manifolds. 

We let (M,g) denote a smooth, compact, connected, n-dimensional Riemannian manifold 
with smooth (possibly empty) boundary dM. If dM = 0, then we assume that the Euler 
characteristic x(M) does not vanish. We study the following system of nonlinear partial 
differential equations: 

u t + V u u = -gradp + z^DivDef u - ADiv(Vd T • Vd) , (1.1a) 

divu(t,ar) = 0, (1.1b) 

dt + V u d = 7 (&d - \(\d\ 2 - l)d) , (1.1c) 



e 2 

u = on <9M, d = h on dM g(h, h) = 1 or dM = , (l.ld) 
u(0, x) = uo, d(0, x) = do and do|aM = h if c?M 7^ 0. (1-1°) 

Here it(t, a;) and d(t, x) are time-dependent vector fields on M, V denotes the Levi-Civita 
covariant derivative associated to the Riemannian metric g, Def u = ^(Vii + Vu T ) denotes 
the (rate of) deformation tensor, the superscript (-) T denotes the transpose, A denotes the 
rough Laplacian of g defined in (2.4), and u, A, 7 are positive constants. In the case that M is 



flat, an open subset of Euclidean space for instance, then V is the componentwise gradient, 
and A is the componentwise Laplacian given in coordinates x l by A = Yl^=i d$dx* ' r ^^ ae 
system of equations (1.1) is the simplified Ericksen-Leslie model |], [?], [| of nematic liquid 



crystals first introduced by F.H. Lin in [pi and later analyzed by F.H. Lin and C. Liu in 
||ll| , O]. Beautiful numerical simulations can be found in fl6| |. 

This system couples the Navier-Stokes (NS) equations with the Ginzburg-Landau (GL) 
penalization of the harmonic map heat flow. The vector field u(t, x) is the velocity field of 
the fluid, while d(t, x) is the penalized (Ginzburg-Landau) approximation to the unit-length 
director field, representing the orientation parameter of the nematic liquid-crystal. 

The parameter e > is the penalization parameter, v denotes the kinematic viscosity of 
the fluid, A is an elastic constant, and 7 is the relaxation-time parameter. 

The coupling term Div(Vd T • Vd) preserves the regularity of the Navier-Stokes equations: 
when velocity u = 0, Ql.laj ) becomes Div(Vd T ■ Vd) = — gradp, so that when dt = 0, d 
is a solution of the static portion of ( |l.ld ) together with the constraint Div(Vd T • Vd) = 
— gradp. This constraint pushes the gradient flow towards a "regular" (H s , s sufficiently 
large) stationary solution. Even though the coupling term has two derivatives, analytically, 
it is essentially identical to the advection term V u u. 

Results. We begin by extending the simplified Ericksen-Leslie model of F.H. Lin to a 
compact Riemannian manifold with boundary. This extension introduces a new curvature 
term in the basic energy laws, and provides a covariant (coordinate-independent) description 
of the liquid crystal dynamics. Motion on the sphere is an important application. 

In Section ^, we prove that the system of equations ( |1.1| ) actually arises from a simple 
variational principle (the system was originally derived using balance laws) . The variational 
principle is the key to our analysis, for it gives the correct scaling; namely, it shows that when 
d is taken to have one derivative greater regularity than u, the liquid crystal system behaves 
as if it were parabolic, a fact which was not previously known (see p|| ). In fact, according 



to [14 1, because the interaction term Div(Vd • Vd) formally has as many derivatives as the 
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diffusion term, the standard Galerkin procedure for obtaining local solutions had failed in 
prior attempts. 

In Section ||, we give a very simple proof of local well-posedness of the system (1.1) on 



(M, g) (in Theorems [T] and ^) using the contraction mapping theorem; this significantly 
simplifies the clever, but lengthy, modified Galerkin procedure employed by Lin and Liu in 
]Tl| . Moreover, the proof does not require use of either the maximum principle or higher- 
order energy laws. 

In Section |3], we show (in Propositions |l|, ^, and ^) that on two-dimensional Riemannian 
manifolds with smooth boundary (possibly empty), there exists an absorbing set for u in 
H 1 and din if 2 . 

In Section |j| we prove the global well-posedness of the system ([Ll]), as well as the 
existence of absorbing sets for u in H s and d in H s+ , and hence of a compact global 
attractor when the dimension is n = 2 (see Theorems [3| and |j) . When dM = brute- force 
energy estimates may be computed, but when dM is not empty, we use the Ladyzhenskaya 
method to obtain the uniform bounds. We remark that the existence of global attractors 
for this system was not previously known. We also remark that since the Navier-Stokes 



equations are a subsystem of (1.1), one does not expect to be able to prove results in 



dimension three which do not already exist for the Navier-Stokes equations; namely, the 



problem of unique classical solutions remains open, while weak solutions exist [11] 



Finally, in Section |g, we introduce the Lagrangian averaged liquid crystal equations (p.6| 



This system is based on the Lagrangian averaged Navier-Stokes equations (see [17| and ref- 
erences therein) , and is derived by averaging the Navier-Stokes flow over small spatial scales 
which are smaller than some positive small number a. We show that this averaged system 
retains the structure of the original system derived by Lin in the form of averaged energy 
laws, but has the advantage of being globally well-posed on three-dimensional domains (see 
Theorem g). The averaged energy law shows that when both the fluid flow is averaged 
together with the director field, both u and d scale similarly, and d is not required to have 
one-derivative greater regularity. Of course, physically, it seems much more natural to us 
to average the fluid flow, since the molecular orientation is already an averaged quantity. 
We believe that the averaged liquid crystal system will be the ideal model for numerical 
computation. 

Some Notation and Interpolation Inequalities. We shall use the notation H S (TM) to 
denote the £P-class vector fields on the manifold M. The H S {TM) inner-product is given, 
in any local chart, by 

s 

(u,v) s = *£(D a u,D a v), 

H=o 

where 

(u,v) = g(x)(u(x),v(x))n(x) 

J M 

denotes the I? inner-product, a = (ax, ...,a n ) is a multi-index, and 

\a\ = ax + ■ • • + a n , D a = d* 1 ■ ■ ■ d% n . 
We shall denote the H S (TM) norm by 



STEVE SHKOLLER 



(•, •) = (•, -)o, and I • I = I • |o- We set Hq(TM) to consist of those vector fields in H (TM) 
which have zero trace on dM. Similarly, vectors in H\(TM) have trace h on dM. We let 



i 



H s h (TM) = H S {TM) n H s ~2(TdM), g(x)(h(x), h(x)) = 1 Vx G dM s > 1 
denote the space of H s vector fields on M which have (H 1 ) trace h on dM and where 
h G H s ~^(TdM). 

For each x G M, we let Bf = {u£ T X M | ff(x)(t>, v) < 5}, and set B s = \J X&M B S X . We let 
H S (M,B S ) denote the iP-class maps from M into -B . 
We have the product rule 



D a (fg)= T c a ,JDPf) (D a -Pg 



\P\<\a\ 
a-(3>0 

For any integer s > 0, we set 

D s u = {D a u : |o| = a}, \\D s u\\ LP = ^ \\D a u\\ LP . 

\a\=s 

We define the spaces 

V = {u£ C°°(TM) | divu = 0, g(u,n) = on 

W={«£ C^(TM) | divu = 0}, 

and through-out the paper, we shall use W s and V s denote the closure in H s of V and W, 
respectively. It follows that 

V s = {ue H S (TM) | divu = 0, g(u,n)| 9M = 0}, 
W s = {ue H S (TM) n Hq(TM) I divu = 0}. 

In section ||[ we shall give an equivalent definition of W s using powers of the Stokes operator. 

We shall need some standard interpolation inequalities, which follow from the Gagliardo- 
Nirenberg inequalities [pc|l , 1 24 ] : 

Suppose 

1 i ( 1 m\ 1 

- = - + a[ +(l-o)- 

p n \r n J q 

where i/m < a < I (if m — i — n/r is an integer > 1, only a < 1 is allowed). Then for 
f :M -» TM, 

idVIlp < ci^/ii, • i/iir (i.2) 

In what follows, we shall use C as a generic constant. Some specific cases in two dimensions 
(n = 2) that we shall need are as follows: 

\v\l°° < C\D 2 v\][i\v\ l li (1.3) 
\v\ l a < C\Dv^\v^ (1.4) 
< C\v^ ltrn \B m v^[T ■ (1-5) 



Equation (O) is often called the Agmon inequality, while (|1.4f )- (ll.5D are often referred 
to as the Ladyzhenskaya inequalities. 
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We shall use div for the divergence operator on vector fields, and Div for the divergence 
operator on sections of T*M <8> TM. 

2. The Variational Principle 
In this section, we shall explain how the system of equations (|0|) arise from a simple 



variational principle, ( 1.1a ) being the first variation of the action with respect to the La- 
grangian flow variable, and ( 1.1c ) being the I? gradient flow of the first variation of the 



action with respect to the director field. It was not previously known that ( |1.1|) can be 
obtained from a variational principle; rather, balance arguments were invoked to derive the 
model. 

We let r](t, x) denote the Lagrangian flow variable, a solution of the differential equation 

dtf](t, x) = u(t, rj(t, x), r/(0, x) = x. 

For / = [0,T], and each t € I, u € C°(I,W S ), s > (n/2) + 2, the map rj(t, ■) : M ->• M 
is an H s volume-preserving diffeomorphism with H s inverse, and restricts to the identity 
map on the boundary dM. We shall denote this set of maps by T> S ^ D . It is a fact, that for 
s > (n/2) + 1, the set T> s D is a C°° (weak) Riemannian manifold (see || and pi] ). 
We define the action function S : V S ^ D x H S+1 (TM) n Hq (TM) -> R by 

S(v,d) = l [ [ {g(v(x))(u(t, V (t,x)),u(t,r,(t,x))) 

1 J I JM 

-Xg(rj(x)) (V [d(t, V (t, x))} , V [d(t, V (t, x))} ) + 2F(d)} y, dt, (2.1) 



where F(d) = ~ l ) ■ Notice that 



where 



f(d) = graded), 
f{d) = 1 (| d |2 _ !) d 



is the (GL) nonlinearity in ( l.lc[ ). The first term on the right-hand-side of ( |2.1| ) is the 



kinetic energy of the fluid, the second term is the elastic energy of the polymers, and the 
third term is the unit-length constraint on the director field d. As a consequence of the 
right-invariance of S with respect to the lifted action of T> s ^ Dl we may compute the kinetic 
energy of the fluid as well as the elastic energy along the particle trajectory rj(t,x). The 
interaction, or coupling, between the velocity u and the director d comes precisely from the 
elastic energy being computed along the Lagrangian flow n(t, x). 

The elastic energy (1/2) f M \Vd\ 2 y is a simplified form of the Oseen-Frank energy, given 
upto the null-Lagrangian by 

/ [ki\ div d\ 2 + K 2 \d x curl d| 2 + K 3 |d • curl d| 2 ] /x (2.2) 

Jm 

The terms in the integrand represent, respectively, the energy due to splay, bending, and 



twisting of the polymers in the nematic liquid crystal. When k = K\ = k 2 = K3, then (2.2) 
reduces to k j M |V<i| 2 ^. We see that in the Eulerian frame, for a director field which is 
exactly taking values in the unit sphere, the energy is given by 



Energy = \ [ (\u(x)\ 2 + A|Vd| 2 ) n. 
1 JM 
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The penalized form of this energy is then 

E=\f {\u(x)\ 2 + \\Vd\ 2 + 2\F(d))v, (2.3) 
where we suppress the explicit dependence on the small parameter e > 0. 



The action (2.1) is the right-translated time- integral of the energy function ( |2.3| ). The 
penalization was motivated by the study of harmonic maps of simply-connected domains VL 
into spheres (see Q); in particular, the space HKQ,^ 1 ) = when |degree(/i)| > 1 so that 
only infinite energy minimizers exist. As a fix for this problem, the penalization method is 
invoked, which enlarges the space of potential minimizers to f^(f2,R 2 ) (which is obviously 
not empty) and simultaneously imposes the unit-length constraint. 

To compute the first variation of S with respect to rj, we let e i— ► <p e be a smooth curve 
in T> s ^ D such that <fro = e, and (d/d£)\ £= o4> £ = w € W s . We let D/de denote the covariant 
derivative along the curve 4> e . Let rf = rj o (f) e so that 

rp = i], and (d/de)\ e =Qrj e = w or]. 

Then (setting A = 1 for the moment), 

d 

de 



S(ri £ ,d)= / { ff (r / (t,x))(( J D/&) e=0 5 t r ? £ (t,x),5 t r ? (t,x)) 

e=0 J I JM 



+ g(rj{t, x)) {{D/de) £=0 V {d{t, </> e (ri(t, x))) , V (d(t, v (t, x)))} dx dt, 

where V is computed with respect to the moving Lagrangian coordinate y = r](t,x), and 
where we have used dx to denote the Riemannian volume- form fi. We use D\ and D2 to 
denote the Frechet derivatives of S with respect to r\ and d, respectively. Integrating by 
parts, and using the fact that dti] = u o r\ and that rj has Jacobian determinant equal to 
one, we see that 

{D 1 S(r j ,d),6r l )= [ [ g(r](t,x)) (-{(D/dt)d t 7](t,x),w(t,r](t,x)) dxdt 

J I JM 

+ 11 9(y)((D/de)\ E=0 Vdo(l> E (y).T ( j> E (y),Vd(t,y))dydt 

J I JM 

9(y) ((-ut(t, y) - V u tt(t, y) - gradp(t, y),w(t, y)) dy dt 

1 JM 



[ [ {g(y)(V w (Vd),Vd)+g(y)(X7d(t,y)-Vw,X7d)}dydt 

J I JM 

9{y) ((-ut(t, y) - V u tt(t, y) - gradp(t, y),w(t, y)) dy dt 



+ 



+ / / 9 (y)(-Div(Vd T -Vd),w)dydt, 

J I JM 

where the last equality follows from the fact that (V lu (Vd), VcZ) = 0, since div w; = 0. Thus, 
since w is an arbitrary variation of rj, we arrive at the Euler-Lagrange equation 

ut + V u m = — gradp — Div(Vd T • Vd). 
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The viscosity (diffusion) term follows from the Ito formula by allowing r](t, x) to be a 
stochastic process, and replacing deterministic time derivatives with stochastic backward- 
in-time mean derivatives (see ||). Thus (|l-la| ) follows as the first variation of the action 
function S with respect to r\. Equation ( |l.lb| ) follows immediately from the fact that rj is 
volume-preserving . 

Letting d £ = d + edd, a much simpler computation verifies that 



(D 2 S( V ,d),5d) = j- 



S( V , d e ) =!! g{y) (Ad - f(d),Sd) dy dt, 



£=0 JUM 



where 

Ad = V*V (2.4) 

is the rough Laplacian and V* is the L 2 formal adjoint of the covariant derivative V. Hence, 
equation (1.1c) is simply the I? gradient flow of d \— > S(rj,d) given by 

j t (d(t, r,(t, x)) = D 2 S( V , d)=Ad- f(d). 

We remark that 

Div(V<F • Vd) = Vd T ■ Ad + g{R{e i: -)d, V ei d), (2.5) 

where R is the Riemannian curvature tensor which is defined for vector fields X, Y, Z on M 
by 

R(X, Y)Z = VxVyZ - VyVxZ + V [x ,y]Z, 

and where {ej} is any local orthonormal frame. The curvature term in equation ( p.5[ ) will 
play an important role in the energy behavior of the system. 



3. Local Well-posedness 
Let P denote the Leray orthogonal projection from L 2 (TM) onto W°, and let 

A = -PDivDef 

denote the Stokes operator, an unbounded, positive, self-adjoint operator on W°, with 
domain D(A) = H 2 (TM) CiW 1 . As usual, we set 

W s = D(Al), s>0. 

This is a Hilbert space with inner-product (A^u, A^v) for u, v G D{A^). The norm \A*u\ 
is equivalent to the H s norm. 

We first prove the local well-posedness of classical solutions. 

Theorem 1. For s > § + 1, and u G W s , d G ^ +1 (TM), there exists T > depending 
only on the data and M , such that 

u G C°([0,T], W s ) , deC° ([0,T],H s h +1 {TM)) 

are solutions to the system of equations ( [OU . 
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Proof. It will be convenient to recast the equations ( |l.lcj ) and (Lie) so that the solution has 
zero trace on dM. For any boundary data h € H s+ z (TdM), we may choose ip £ H s+ (TM) 
such that trace^) = h. Let 

d = d — tp, so that d\oM = 0. 
We rewrite the system ( |l.l|) as an evolution equation in 

X s = W S ®H S Q +1 (TM) : 



u t + vAu + PV u u = -PXBw(V[d + ^f ■ V[d + ^]), (3.1a) 
dt + V u d = 7 (Ad - /(J) + Av) , (3.1b) 
it = on dM, d = on dM or dM = , (3.1c) 
u(0, x) = no, d(0, x) = do(x) = do(x) + ip(x) , (3. Id) 

where 

/(rf) = l(|d + Vf-l) + 

We define the vector 

x = (u,d) € X s ; 

since i2" s_1 -class vector fields form a Schauder ring for s > \ + 1, we may define the maps 

(3.2) 



2 

1 : X s -> y 8 " 1 , 



= -P (V u u + Div(V[d~ + • V[d~ + , 
and 

4> 2 ■ X s -» H S (TM), 
Mx) = -Vj--yf(d) + 7 AV. (3 ' 3) 

Thus, the vector 

* = (01,02) : X s -> V^" 1 x H S (TM). 
We are using the fact that the projector P maps H s ~ l to itself. To see this, we write 

PDiv(Vd T • Vd) = Div(Vd T • Vd) - gradg, 

where q solves the Neumann problem 

Aq = divDiv(Vd T • Vd), 
g(grad q,n) = g (Div( Vd T • Vd) , n) on dM. 

Since div Div( Vd T - Vd) is in H S ~ 2 (M) and Trace Div(Vd T - Vd) is in H s ~^(dM), by elliptic 
regularity (/ is in H s (M) so that grad (/ is in H S - X {TM), as desired. One sees that PV u u is 
also in H s (TM) by a similar argument. 
Next, we define the semigroup 

- -tvA Q 

*(«>= „ e nA 
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We can now express the system ( |3.1[ ) as the integral equation 

Xt (t,-) = S(t)x Q - [ S{t- s)$(x{s))ds = tfx(V). (3.4) 
Jo 

Since e~ tvA : W s -> W s and e^ A : H° +1 (TM) -> H s Q +l {TM) are strongly continuous 
semigroups, it follows that 

£(t) : X s — » X s is a strongly continuous semigroup for t > 0, (3-5) 

and that for i > 0, S'(t) : V s-1 x H S (TM) — > X s ; furthermore, we have the usual estimate 
(see, for example, f24fl ) 



(3.6) 



(3.7) 



\l{V s - 1 xH s (TM),X>>) <Ct 2, (0,1]. 

Using the fact that for s > (n/2) + 1, P : H S (TM) n H s ~^{TdM) -> V s is a bounded 
projection, we obtain that 

$ : X s -> y s_1 x H S (TM) is a locally Lipschitz map; 

namely, 

|| 0i (it, J) - 0(u,e)|| s _i < Ci (||u - v|| a , || J- e|| s+ i 
||(/> 2 (tt,d) - 0(f,e)|| s < C 2 M|ii - v\\ a , \\d- e|| a+ i 

where Ci and C2 depend on ||u|| s , |M| S , ||<i|| s +i, ||e|| s+ i, and HV'lls- 
Fix a > and set 

Z = {x€C([0,T],X s ) I x(0) = (n , d ), \\x(t, ■) - x(0)\\ x * < «}• 

We want to choose T sufficiently small so that \l/ : Z — > Z is a contraction. By ( |3.5| ), we 
can choose Ti so that 

||5(t)x -a;o||jC« < «/2 Vt€[0,Ti]. 
If x € Z, then by (^) we have a bound 

||*(x(s))||v-ixH'(TW) < K\ f °r b € [0,Ti]. 
Using (|3lf ), we have that 



S(t-s)$(x(s))ds 

hence, for i € [0, T2], and with x = (u,d), 

\\y(u(t),d(t))-y(v(t),e(t))\\ X s 



S(t-s) [$(u(s),d(s))-$(v(s),e(s)) 
< CtiK sup \\{u(s),d(s)) - (v(s),e(s))\\ X s- 



ds 



Choosing T < T 2 small enough so that CT2K < 1, we see that by the contraction mapping 
theorem, ^ has a unique fixed point in Z, and this proves the theorem. □ 

Using the contraction mapping theorem, we can also establish the local well-posedness 
for a weaker class of solutions. 
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Theorem 2. Suppose 2 < dim(M) < 5 and set sq = § + |- For s G (so, 2) and no G VK S , 
do £ H^ l {TM), there exists T > depending only on the data and M , such that 

n G C° ([0, T], W s ) , deC° ([0,T],H s h +1 (TM)) 

are solutions to the system of equations \1. 

Proof. We keep the same notation as in the proof of Theorem p]. For s G (so,2), we have 
that 

\\S(t)\\L(y°xHi<TW) < Ct~\ 7 e (0, l), t g (0, l]. 
Thus, it suffices to prove that for s G (s ,2), the map <E> : X s — > 1/° x H 1 (TM) is locally 
Lipschitz. Using Lemma 5.3 [[24], Chapter 17], we have that for s G (so, 2) (/,<?) i— > / g : 
H s x H s —> H 1 . It follows that u i-> n <g> n : i? s -» iT 1 , d ^ (Vd T • Vd) : # s+1 -» iT 1 , and 
(n,d) i — > V u d : H s x iT s+1 -» H 1 . The fact that d i-> /(d) : iT s+1 -» H 1 follows because 
.£P +1 forms a Schauder ring. Hence, $ is indeed locally Lipschitz, and the remainder of the 
proof is identical to the one for Theorem |l|. □ 

4. Basic Energy Laws on Riemannian manifolds 



In this section, we show that the system ( |1.1|) admits the following energy law: 
d 



E = - (i/|Defn| z + A7|Ad-/(d)|^) - A Trace(i?(-, n)d, V.d), (4.1) 



where E is given by and ej denotes a local orthonormal frame. When M has zero 

curvature, then E is a Lyapunov function for the system (LI), with the property that 

E{u(t),d(t)) < E(u 0) d ), Vt > 0, 

and if E(u(ti),d(t\)) = E(u(t 2 ), d(t 2 )) for t\ < t 2 , then (u(t),d(t)) = (n*,d*) are equilib- 
rium solutions. Even, when the curvature R ^ 0, the energy remains uniformly bounded. 

This bound, in turn, then yields an a priori uniform bound for the pair (n, d) which shows 
that all solutions eventually enter an absorbing ball in W 1 x H 2 (TM). 

The following two lemmas are standard: 

Lemma 1. If dM / 0, then for d G H\, 

|d(t,-)k« < li vt > 0. 

Proof. We compute the pointwise inner-product of ( [4,2b| ) with d, and use the fact that 
#(Ad(x),d(x)) = (l/2)A(g(d(x),d(x))) - o(Vd,Vd). Hence, Vx G M, we obtain 

i^o(d, d) + i«(grad[<7(d, d)], n) - ^A(g(d, d)) + o(Vd, Vd) = ~\g(d, df - g(d, d)]. 

Now suppose that max tjX g(d(t, x), d(t, x)) occurs at (to,xo), an element of the parabolic 
interior; then 

^g(d(t ,x ), d(t ,x )) = 0, Vd(t ,x ) = 0, grad[g(d(i , x ) , d(t , x ))\ = 0, 
and Hess g(d(to, xo), d(to, xq)) < 0. This implies that 

-^A(g(d{t , x ), d(t , x ))) > 0, 
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but 

--^[g(d(t ,x ),d(t ,x )) 2 - g(d(t ,x ),d(t ,x ))} < 0, 
which is a contradiction. □ 

Lemma 2. If dM = 0, then for d € H 2 (M, B 8 ), 

\d(t,-)\ L °° <S, Vt > 0. 
Proof. This again follows from the maximum principle above. □ 



Proposition 1. The energy law l \4-V holds, and there exists an absorbing set for (u,d) £ 
W° x Hl(TM) if dM f and for (u, d) € W° x H X (M, B s ) if dM = 0. 

Proof. Using the formula ( |2.5| ), we rewrite (1.1a) and ( |l.lcj ) as 

u t + V u u = -gradp + i^DivDef u - XVd T ■ Ad - g(R(ei, -)d,V ei d) , (4.2a) 



1 



dt + V u d = 7 Ad - ^(|d| 2 - l)d) , (4.2b) 



c 



ld_ 

2dt 



where e\ is any local orthonormal frame. Adding the L 2 inner-product of ( [4.2a ) with u to 
the L inner- product of flPbD , we obtain the basic energy law 

u\ 2 + A|Vd| 2 + 2A / F(d(x))fi) 

= - (y\T)eiu\ 2 + \j\Ad- f(d)\ 2 ^j — A Trace(i?(-, u)d, V.d). (4.3) 

In the case of a flat manifold, such as a bounded domain in W\ R = 0, and ( [4.3j ) reduces 
to the basic energy law (1.8) in [11 1. 
From Lemmas |l| and ^ we have that 

\d(t,-)\ L oo <C, t>0. (4.4) 

It follows that 

Trace(i?(-,u)d, V.d) < C\R\ L °°\u\\d\^\Ad\^ 

< Ce\Ad\ 2 + ^ (\M\\R\ L °o\u\)% 



< Ce\Ad\ 2 + e\u\ 2 + -j (|M||«| LO o) 4 

< Ce|Ad| 2 + Cb 1 (M)e| Def u| 2 + ^ (|Mj|i?j L 



V 1 



where the second and third inequalities follow from Young's inequality ( |5.2| ), and the last 
inequality follows from the Poincare inequality for cq{M) > 0, a positive constant depending 
on M. Taking e > sufficiently small so that 

K = min(co — e, 1 — 2e) > 0, 



the basic energy law (4.3) on a Riemannian manifold yields the following differential in- 
equalities: 
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ld_ 
2~dt 



\d_ 

2dt 



\u\ 2 + \Vd\ 2 + 2 f F(d)fi 

JM 

<-K C |Defu| 2 + |Adj 2 4 
\u\ 2 + \Vd\ 2 + 2 f F{d)n 

JM 

\u\ 2 4 \Vd\ 2 + 2 / F(d)/x 
ijif 



< -If C7 



4 Po 



(4.5a) 



(4.5b) 



where 



Po = C 

Using the classical Gronwall lemma, we obtain 
|u| 2 + |Vd| 2 + 2 f F(d)n 

JM 



(K + C/e - l)\M\ + —\M\ A \R\ L oo 



< 



Kl 2 + | Vrfol 2 + 2 / F(d V 



+ PD (l- e -*<*) 



Thus, 



lim sup 



t— >oo 



K*)| 2 + |Vd(t)| 2 + 2 / F(d(t))/i 



< Po 



lim sup [|ii(t)| 2 + \Vd(t)\ 2 ] <p + 2\M\ 



(4.6a) 
(4.6b) 
□ 



When R = 0, we do not need to rely on the maximum principle to establish Proposition 
H or to establish the existence of an L°° absorbing set for ( 4.2b| ). 

Lemma 3. If M = and R = 0, then for do £ H 2 (TM), there exists > and some 
t* > independent of do such that 

\d(t,-)\ L °o < Poo Vt > t*. 

Proof. When R = 0, from the energy law ( |4.1| ), we see that there exists and L 2 absorbing 
set, so that for some t > t* all bounded subsets of L 2 (TM) will enter the L 2 ball of radius 
Po- 

For p > 2, we take the pointwise inner-product of ( 4.2b| ) with p\d\ p ~ 2 d and integrate over 
M to obtain the differential inequality 

±\d\ p LP = -p[ \vd\ 2 \dr 2 p-p( P -2) [ \dr 2 \v\d\\ 2 p 

at JM JM 



+ ji [ \d\ P LP ~ \d\ p L ~ P 2 



<- P (p-2) [ \ d r 2 \v\d\\ 2 p+\\d\ p LP . 

JM e 



(4.7) 
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Using the interpolation inequality (see [15| for details and further applications) 



p-2 

\d\i P <c p \d\ 2 (Jjv\d\ p *\ 2 ^ p =c p \d\ 2 (^jjdr 2 \v\d\\ 2 i?j p 

we see that 



P-2 



-pip -2) / \d\v~ 2 \V\d\\ 2 LL < 



1 \ W^l 4p(p - 2) & 



Using Bernoulli's trick in the differential inequality Q4.7| ), we get a uniform bound for 
\d(t, -)\lp which is independent of p (even if the constant C p tends to infinity), and thus we 
may pass to the limit as p — > oo. □ 

Using Lemma [| we immediately have 

Proposition 2. If R = and dM = $, we have the energy law 

^|Defu| 2 + A7|Ad-/(d)| 2 ) , 
and there exists an absorbing set for (u,d) € W° x H 1 (TM). 



dt 



Proposition 3. For dim(M) = 2, there exists an absorbing set for (u,d) 6 W 1 x H?(TM) 
if dM / 0, for (u,d) £ W 1 x H 2 (M,B 5 ) if dM = and ^ 0, and for (u,d) £ W 1 x 
H 2 (TM) if dM = and R = 0. 



Proof. It follows from (|4.5a|) that 

rt+r 



KC 



SO 

lim sup 

t^oo Jt 

Therefore, 

rt+r 



Def u(s)\ z + |Ad(s)r + 2 / F{d{s))n 

'm 



ds 



< rp + 



|u| 2 + |Vd| 2 + 2 f F(d))n 
Jm 



t+r 



Defu(s)\ z + \Ad(s)\ 2 + 2 / F(d(a))/x 

'm 



Vr > 0, 
ds < (r + l)p - 



|Defu(s)| z + |Ad(s)| 2 + 2 / F(d(s))(i 



ds is uniformly bounded. 



(4. 



Now let 

,4 2 = |Defu| 2 + |Ad-/(d)| 2 , 5 2 = |VDef u\ 2 + |V(Ad - f(d))\ 2 . 
Using (^8|) and flO]) , we have that 



t+r 



A 2 (s)ds is uniformly bounded. 



(4.9) 



In the case that R = 0, it follows from a similar argument as in (4.4)-(4.8) of |1 .] that for 
some constants c\, 02,03 > 0, 



d A 2 (t) + Cl B 2 {t) <c 2 A\t) + c 3 . 



dt 



(4.10) 
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When R / 0, we find that for C4 > 0, 

j t A 2 {t) + ciB 2 {t) < c 2 A A {t) + c 3 + c 4 Trace( J R(-, Au)d, V.d). 

The last term is bounded by ej-Rl^oo | Div Def u\ 2 + (C / e)\V d\ 2 , so by taking e > sufficiently 
small and adjusting the constants as necessary, we see that flllCD still holds. 

Thus, using (4.6) and appealing to the uniform Gronwall lemma (see, for example, p3[ ), 
we see that 

j4(t) is uniformly bounded in time. 

Because of (|4.4| ), we may extract a uniform bound for | Def u\ 2 + |Ad| 2 . Hence, we have an 
a priori uniform bound for u in the H l topology and for d in the H 2 topology. □ 

5. Global Well-posedness and Global Attractors 

We shall first consider a closed Riemannian manifold such as, for example, the two-sphere 
§ 2 ; for such manifolds, simple brute-force energy estimates work. 

Theorem 3. Forn = 2, s> I, dM = 0, and u E W s , d E H S+1 (M, B 5 ), 

u £ C°([0,oo],W s ), d E C°([0,oc},H s+1 (M,B 5 )) 

are solutions to the system of equations 1\ ). Moreover, there exists a compact global 
attractor for the system in VF 5-1 x H S (M, B 5 ). In the case that R = 0, we can replace 
H S+1 (M,B 5 ) with H s+l {TM). 

Proof. Taking the H s inner-product of ( p.. la ) with u and adding the H s+1 inner-product of 



^ (Ml + \d\ 2 +1 ) < -u\u\ 2 s+1 - 7 |d|? +2 + A|(PDiv(V(F • Vd),u) 



(1.1c) with d, we find that 
ld_ 
Tdt 

+ I (PV u u, u} 3 \ + I (V u d, d) s+1 1 + 7 I (/(d), d) s+1 1 . (5.1) 

We shall estimate each of the nonlinear terms on the right-hand-side of ([Tl]); as we showed 
in the proof of Theorem |l], the projection P acting on the nonlinear terms, maps H s ~ 1 into 
itself continuously, so it suffices to estimate (\7 u u,u) s and (Div(Vd T • Vd),u) s in the third 
and fourth terms. Using Proposition [|, we may interpolate the nonlinear terms between 
\u\i and |u| a +i and \d\2 and |d| s +2, respectively. 
We have that 

(V u u,u) s = ^2{D a (V u u),D a u) =J2 E {caAD P Vu)(D a -Pu),D a u) 

a=s ot=s \p\<e 

a-/3>0 



E E MDNu){D<*-Pu),D a u) 



i \0\<S-1 

a-/3>0 



1 

< C E l^+^l^lD^l^l^ul 

m=0 



where we set = m so that |a — (3\ = s — m, and the last equality follows from the fact 



that (S7 u (D s u), D s u) = 0, since divu = 0. For m = 0, ...,s — 1, we use ( |l.4|) and (|l.5[) to 
estimate 
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\D m+l u\ Li \D s ~ m \ Li \D s u\ < C\D m+1 u\^\D m+2 u\^\D s ~ m u\^\D s ~ m+1 u\^\D s u\ 

3+1 2(a-l) 

<C\u\^ \u\ s £ . 
Using Young's inequality, 



X C 1 

a x b<ea + —b^, a,b>0, < A < 1, (5.2) 



it follows that 

C 

(PV u u,u) s < e\u\ 2 s+l + -\u\{ +1 . 
For the next term, we have that 

(Div(VcF • Vd),u) s = {D a Div(Vd T • Vd),D a 

\a\=s 



^ C J2 J2 I (D (3 Vd)(D^I 3 Vd)(D a u) f i 
, fit Jm 



\<x\=b \/3\< 
| 7 |=s+l7-/3>0 

s+1 



< C ({D m+1 d)(D s - m+2 d), {D s u)), (5.3) 

m=0 

where m = \(3\. In the case that m = 1, s, ( |5.3[ ) is bounded by 

C £ \D m+1 d\l 4 \D s ~ m+1 d\l 4 < C\d\ 2 \d s+2 \u\\ Mjg 

m=l 



<s\d\j +2 + j\d\i\u\]\u\ 8+ i , (5.4) 



where the first inequality follows from repeated use of (L5), and the last inequality follows 
from ab < ea? + (C/e)b 2 , where a, b > 0. One more application of ([)!]) shows that ( |5.4D is 
bounded by 

C 

e|c^ +2 + e\u\ 2 s+1 + ^+jM| s |u|?. 



In the case that m = 0, s + 1, ( |5.3| ) is bounded by 



ClciU+alrfU* |X>*t*|x,4 < C\d\ s+2 \d\ 2 \u\l°\u[^ 



\s+l 

(J 1 2(s-l) 

<e\d\ 2 s+2 + -\d\ 2 2 \u\{ \u\ s+ { 

C 

< e\d\ 2 s+2 + e\u\ 2 s+1 + —IdH-lul!, 

where the first inequality follows from (|ll]), and the last two inequalities follow from ( [5.2| ). 
It follows that 

(PBiv(Vd T ■ Vd),u) s < e\u\ 2 s+1 + e\d\ 2 s+2 + ^\d\ 2s (|«|i + |«|f) . 
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We next compute that 

(V u d,d) s+1 = £ (D a (V u d),D a u) 

\a\=s+l 

<C Yl ({D (3 Vd)(D a ~ l3 u),D a Vd) 

|a|=s+l \/3\<s+i 
a-f3>0 

s 

< C ({D m+1 d)(D s ~ m+1 u),D s+1 d), (5.5) 

m=0 

since for m = s + 1, we have that (V u (D s+1 d), D s+l d) = 0. We estimate the case m = 
first in (|5.5|): 



\Dd\ Li \D s+1 d\ Li \D s+1 u\ < C\d\ 2 \u\ s+1 \d\] +1 \d\* +2 

2s+l 



<C\d\ 2 *° \u\ s+1 \d\s 2s 



<e\d\ 2 s+2 + -( [ \d\ 2 * \u\ s+1 ) 

<e\d\l +2 + e\u\ 2 s+l + ^\d\\ s+ \ 

where the last two inequalities follow from two applications of the Young's inequality. 

For the cases 1 < m < s, is bounded by C\D m+1 d\ L A \D s ~ m+1 u\ L A \d\ s +\, so by Q 
and ( |1.5| ), we find that for m = 1, s, 

\D m+1 d\^\D m+2 d\^\D s - rn+l u\^\D s - m+2 u\^\D s+1 d\ 

2s-2m+3 2?re+2s-3 2?n-l 2s-2m + l 

<c\d\ 2 - |d[ 8+ /- \u\ s+ *° 



((-i 4m-2 \ / 4s-4?ti+2 

where we have used Young's inequality for the last step. Another application of Young's 
inequality yields the estimate 

(V u d, d) s+1 < e\d\ 2 s+2 + e\u\ 2 +1 + 2B ° 2m+5 \d\f^~ 3 \u\ 2 ™-\ m = 1, s. 

g 2s-2m+3 

For the final nonlinear term, we have that 

s+1 m 

(f(d),d) s+1 ^((^d)(I} m ~ n d)p s+1 ~ m d),# s+1 d} 

m=0 n=0 

< \D n d\ L s \D m - n d\ L s \D s+1 ~ m d\ L 4 \D s+1 d\. 

Using the estimate 

\v\ L s <\v\l 2 \v\l, (5.6) 
together with (|1.5| ) and Young's inequality, we have that 

2(s+2) 2(s-2) fl 

{f(d),d) s+1 <C\d\ 2 ' |d| s+ l <s\d\ 2 s+2 + 7 \d\ s + 2 . 
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Letting 



p = C 



\ {\ d \^ + luir- 1 ) + ^\d\i s (\ u \ 1+ \u\d + ^mr 2 



+ 



y ; Irfia.-am+si |2m-l 

/ ^ is — 2 m + 5 I Iz 111 



^ £ 2s-2m+3 



and taking e > sufficiently small so that 

if = min (z; - 4e, 7 - 4e) > 0, 



the basic inequality (5.1) takes the form 
d 



dt 



(\u\ 2 + \d\ 2 s+1 ) <-K (\u\ 2 s+1 + \d\ 2 s+2 )+Cp 
<-K{\u\ 2 + \d\ 2 s+1 )+Cp. 
Letting c\ = CK > 0, the classical Gronwall lemma gives 



so that 



+ |€+i) < (l«ol5 + \d \ 2 s+ i) e~ Clt + Cp(l- e- C1 '), 



limsup(|u|2 + |d|2 +1 ) <Cp. 

t— >oo 



(5.7) 



Thus, since the time interval of existence from Theorem |l] only depends on the initial 
data, the a priori bound ( |5.7| ) together with the continuation property gives the global well- 
posedness result. Moreover, because of the absorbing sets that exist in W s x H s+l {TM) 
by virtue of (|5.7|), we obtain using Theorem 1. 1.1 of [ 23 1 , the global attractor that we 
asserted. □ 

For a Riemannian manifold with boundary, the above (brute force) H s energy estimate 
does not work, because boundary terms arising from integration by parts on the diffusion 
term v Div Def u do not vanish. It is possible, however, to obtain estimates on ut and dt 
which provide the global well-posedness result. 

We have that 

u tt = v Div Def u t — V Ut u — V u u t — gradp 4 — XV df ■ Ad — XV d T ■ Ad t 



Xg(R( ei , -)d t ,V ei d) - Xg{R(e h -)d, V ei d t ) 



and 



dtt = ~V Ut d - V u d t + ^Ad t - 7 grad f(d) ■ d t . 

Since d t = on dM we see that (|ut| 2 + | Vd t | 2 ) = (u t ,uu) — (Adt,d t t)- Standard 
interpolation combined with Young's inequality yields, for constants c\,C2 > 0, 

7^(M 2 + |W 4 | 2 ) =ci (|Defn| 2 + |Ad| 2 + c 2 ) (\u t \ 2 + \V d t \ 2 ) . 
By Proposition 0, for each t, 



u(t, •) G W 1 and d(t, •) G H 2 (TM). 



(5.8) 
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It follows that if u G W 2 , d G H 3 (TM), and h G H^(TdM), then u t (0) G L 2 (TM) and 
d*(0) G Hl{TM) so that 

Uf G L°°((0,oo), W°) and d t G L°°((0, oo), Hq(TM)). (5.9) 

From ( |5.8| ), we claim that \7 u d is in H 5 (TM) for 5 G (0, To see this, note that for e > 

w^w w : H p ^ H e{1+e \ where P=\ + ^C 1 + e ) + 

We set <5 = 0(1 + e), and, for example, set e = | and < ^; then 5 G (0, ^) and p > ||, 
so the claim is established. Using standard elliptic regularity on equation ( |4.2b| ) , we see 
that d G H 2+ \TM), and the H 2+s -norm of d only depends on the initial data and M. 
This shows that Div(Vd T • Vd) is in L 2 so that with ( |5.9|) , we see that u is in W 2 . By 
bootstrapping, we find that d is in H 3 , and the continuation argument shows that the 
unique solution may be continued for all time. If h G C°°(TdM), then both u and d are in 
C°°((0, oo) x M). Proposition |3] together with Theorem 1. 1.1 of |^] proves the existence of 
the global attractor in W° x Hj^TM). Thus, we have the following 

Theorem 4. Suppose thatuo G W 2 and do G H^iTM). Then there exists a unique solution 
u G L°°((0, oo), W 2 ) and d G L°°((0, oo), F|(TM)). 

If h G C°°{TdM), then both u and d are in C°°((0,oo) x M). Furthermore, there exists a 
compact global attractor in W° x H^{TM). 



6. LAGRANGIAN AVERAGED LIQUID CRYSTALS 

As we described in the introduction, the director field d describes locally the averaged 
direction of the constituent molecules; it is thus reasonable, and of practical and computa- 
tional importance, to locally average the Navier-Stokes fluid motion as well. Recently, the 
Lagrangian averaged Navier-Stokes (LANS) equations were introduced as a model for the 
large scale Navier-Stokes fluid motion which averages or filters over the small, computation- 



ally unresolvable spatial scales (see [17] and the references therein). The LANS equations 
are parameterized by a small spatial scale a > - fluid motion at spatial scales smaller than 
a is averaged or filtered-out. There are two types of Lagrangian averaged Navier-Stokes 
equations: the isotropic and the anisotropic versions. We shall begin with the isotropic 
theory, and for simplicity of presentation, we shall assume that M is flat. 
The isotropic LANS equations for the mean velocity u(t, x) are given by 

d t {l - a 2 A)u + V u (l - a 2 A)u - a 2 Vu T ■ Au = - gradp + u(l - a2A)Au (6.1a) 

div u(t,x) = 0, (6.1b) 

u = 0on<9M, (6.1c) 

u(0, x) = uq. (6. Id) 



Equation ( 6.1a ) has an equivalent representation as 

d t u + V u u+U a {u) = -(1 -a 2 A) _1 gradp + i/Au (6.2a) 
U a {u) = a 2 {\ - a 2 Ay 1 Div [Vu • V/ + Vn • Vu - Vu T • Vu] . (6.2b) 
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When dM = the LANS equations take on a particularly familiar "sub-grid-stress" form 
with fl6,2|) becoming 

dtu + V u u + Div r a (u) = — gradp + uAu (6.3a) 
T a (u) = a 2 (l - a 2 A)" 1 [Vu • Vu T + Vu • Vu - Vu T ■ Vu) , (6.3b) 

where r Q representing the sub-grid or "Reynolds stress." 

The remarkable feature of the LANS equations is that, unlike the Reynolds averaged 
Navier-Stokes (RANS) equations or Large Eddy Simulation (LES) models of turbulence, 
no additional dissipation is put into the system. In fact, when v = 0, the LANS equations 
conserve the Hamiltonian structure of the Euler equations with both a modified kinetic 
energy 



E a -- 

2 



f (\u\ 2 + 2a 2 |Defu| 2 ) fi (6.4) 
Jm 



and helicity 

H a = wAdw, w = (l-a 2 A)u b , u*=g(u,-), (6.5) 
Jm 

being conserved. 

This is easiest to see from equation ( |6.3| ), where the only term that is added (to the 
NS equations) is Divr a (n); it is precisely this term which averages the small scales, and 
this is accomplished by the use of nonlinear dispersion as opposed to dissipation. A simple 
computation, which requires taking the I? inner-product of the LANS equations with u 
when v is set to zero, shows that (6.4) is conserved. Why is it so important not to over- 
dissipate the NS equations? The answer is twofold: first, the addition of artificial dissipation 
obviously and spuriously removes crucial small-scale features, and second, artificial viscosity, 
which is present in RANS or LES models, suppresses intermittency, a fundamental feature 
of fluid turbulence. 

Mathematically, for all a > 0, the three-dimensional LANS equations are globally well- 



posed (see [18]), yet when the averaging parameter a is taken sufficiently small, computa- 
tional simulations of LANS are statistically indistinguishable from the simulations of the 
NS equations. Furthermore, the LANS equations provide a tremendous computational sav- 
ings as shown in simulations of both forced and decaying turbulence ( [0j , |]l9| ) . Finally, the 
LANS equations arise from a variational principle in the same fashion as the NS equations. 
We shall therefore base our development of the averaged liquid crystal equations on the 
LANS model, and introduce the following system of equations: 

u t - uAu + V u u - U a (u) = -(1 - ^A)" 1 [gradp + Div(V(f T • Vd)] , (6.6a) 

dwu{t,x) = 0, (6.6b) 

dt + V u d = i [Ad - \(\d\ 2 - l)d) , (6.6c) 



e 2 

u = 0on dM, d = h on dM g(h, h) = 1 or dM = , (6.6d) 

u(0, x) = uq, d(0, x) = do and cfolaM = h if dM / , (6.6e) 
where A denotes the componentwise Laplacian. 
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Averaged Variational Principle. Following the notation of Section |2[ we define the 
averaged action function S a : V s D x H s+l (TM) n Hq(TM) -> R by 



= ^ ^ Li {9 

+ 2a 2 g(i](t, x)) (Def u(t, 77(4, x)), Def u(t, rj(t, x))) 



2 



+A 5 (r ? (t, x)) (V [d(t, 7?(t, x))] , V [d(t, rj(t, x))]) + 2F(d)} p dt, (6.7) 

where we suppress the explicit dependence of u and d on a and e. 

Again, we see that ( 6.6a ) arises as the first variation of the action function S a with 
respect to 77, and the remaining equations are identical to the original system (|1 . 1|) . Note, 
however, that now u is the mean velocity, and it is the mean flow rj which is transporting 
the director field d. 

Averaged Energy Law. For simplicity, we shall present the formulation in the case that 
dM = 0, although the more general case follows in the same fashion as we presented above. 
Following the notation of Section 0, we have the following basic averaged energy law: 



1 d 
2dt 



(\u\ 2 + a 2 \Vu\ 2 + X\Vd\ 2 + 2 j 



F(d)fi 



< -v (|Vn| 2 + a 2 \Au\ 2 ) - 7 A|Ad - f(d)\ 2 . (6.8) 

From Lemma ||| and ( |6.8[ ), it follows that there exists i > and some po > which are 
independent of the initial data such that 

\u(t,-)\l + \d\l<p Vt>i. 

We see that the averaged energy law is, in some sense, more natural than the standard basic 
energy law (4.1) since the director field d is no longer constrained to have one derivative 
greater regularity than the velocity of the fluid u: both u and d now scale similarly. 

Because of the a priori uniform bound of u(t, ■) in W 1 , it is very easy to obtain an a priori 
bound for u(t, •) in W 2 and d € H 2 {TM) when the dim(M) = 3. We simply compute the 
sum 

= ((1 - a 2 A) (ra) , (1 - a 2 )u) + (A(H), Ad), 

and use Lemma |3| Using similar estimates as above, we obtain an a priori energy estimate, 
in fact an absorbing set, in W 2 x H 2 (TM), and by bootstrapping, we may easily obtain 
higher-order a priori estimates. Local well-posedness follows again from the contraction 
mapping argument that we gave in Theorem |l|, so we have the following 

Theorem 5. For n = 2, 3, s > I, dM = 0, R = and u G W s , d G H S+1 (TM), 

u G C°([0, 00], W s ), d G C°([0, 00], H S+1 (TM) 

are solutions to the system of equations ( \1. Moreover, there exists a compact global 
attractor for the system ( \T7\ ) in W s ~ l x H S (TM). 

It is not difficult to generalize this Theorem to manifolds with boundary following the 
method in 18]. 
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7. Concluding remarks 

Gradient flow versus damping. We considered the 1? gradient flow of the variation 
of the action function S(rj,d) with respect to d in the director field equation ( |l.lc| ). In 
the liquid crystal literature, however, it is common to see a damped second-order equation 
for the director field (see 0j and references therein) , which in the context of our simplified 
system would mean replacing d t with (5\d tt + fodt for some constants 0i and Of course, 
both types of equations have the identical stationary solutions, but in terms of stability, 



L. Simon's result [22] guarantees that the damping term takes over. As far as parabolic 
estimates are concerned, it is easy to treat either type of equation, but we feel it is more 
natural to take the path of steepest descent in relaxing the orientation towards its preferred 
configuration. 

Lie advection versus parallel transport. This remark concerns the coupling term V u d 
in equation ( |l . l<j ) . This term arises by considering the time derivative of (d o rj)(t,x) := 
d(t, rj(t, x)), where for each t, rj(t, •) is a volume-preserving diffeomorphism in the topological 
group V s D . In group-theoretic language, this suggests that the action of V s ^ D on the vector 
space of director fields is on the right. The natural action of T> s ^ D on the vector space of 
director fields, however, is on the left, or by push-forward: instead of dor/, the natural 
action is r]*d := Dr\ • d o if. Taking the time derivative of [r]*d](t,x) gives C u d, the Lie 
derivative of d in the direction u. The Lie derivative C u d = V u d — VdU, and this is the 
actual term which is present in the Ericksen-Leslie model. A nontrivial extension of our 
analysis is required to analyze the system (LI) with V u d replaced by C u d, and this shall 
be the focus of a future article. 

Other fluids models. Using our methodology, it is quite easy to study a number of other 
fluids models. For example, by replacing the Oseen-Frank energy with the Landau-Lifshitz 
free energy 1/2 J M Ag(V M, V M)/j, where M is the direction of magnetization in a cubic 
ferromagnet, we can obtain an almost identical system of PDEs. Similarly, if we replace 
the vector d in our action function S with a scalar field <fi, and replace the L 2 gradient 
flow in equation ( |1 . lcj ) with H -1 gradient flow, we obtain a model of two-phase flow whose 
interface moves via motion by mean curvature (see fl5|]). This model consists of a coupled 
Navier-Stokes Cahn-Hilliard system, where the interface is governed by surface tension. 

The defect law in the limit as e — > 0. We considered the GL penalization of the Oseen- 
Frank energy law so as to obtain finite-energy minimizers, but we have yet to consider the 
limit of our solutions as e — > 0. It remains an open problem to characterize the dynamical 
law of the GL vortices when coupled to the Navier-Stokes motion. Following the pioneering 



work in ||] and [13] on the dynamical law of the GL vortices, we expect that the location 
of the jth vortex, a,j, will solve the distributional equation 

j- t \ sw 

OtCLj + aiv(ajU) 



s aj 



where u simultaneously solves the Navier-Stokes equations, and W = Yltyj 1°S I 

X <i X j | IS 

the renormalized energy. We expect that a rigorous defect law will be much easier to obtain 
when u is instead a solution of the LANS equations, because in that case, u is uniformly in 
H 1 with respect to the penalization parameter e > 0. 
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